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We present results from numerical simulations of the 2+ld SU{2)®SU{2) Nambu-Jona-Lasinio 
model with Nj = 4 fermion flavors at zero and nonzero temperature T. At zero temperature, 
critical exponents extracted from the scaling of the order parameter and fermion mass are found 
to be consistent with next-to-leading order predictions of the 1/iV/ expansion. At nonzero 
temperature we observe fermion mass generation despite the lack of chiral symmetry breaking, 
which is forbidden by the Colemann-Mermin- Wagner theorem for all T > 0. We study the 
effects of lattice discretisation and finite volume on the dynamically generated fermion mass. 
By studying the lattice dispersion relation we also show that in the hot phase there is no 
significant temperature induced modification to the speed of light. Studies of the equation 
of state are made by measuring the pressure as a function of temperature and comparison is 
made with large- A'^^ predictions. 
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1 Introduction 



The nature of QCD and related field theories at finite temperature remains a highly 
active area of both analytic and numerical research. By studying the high temperature 
transition in strongly interacting theories we gain insight not only into the hot chi- 
rally symmetric phase, but also the cold phase, in which spontaneous chiral symmetry 
breaking leads to a dynamically generated fermion mass. 

In 2 + 1 dimensional field theories the nature of the transition can be particularly 
interesting. It has been shown 0, |^ that in the three-dimensional Gross-Neveu model 
with a U{1) chiral symmetry at nonzero temperature, there is a regime at tempera- 
tures above the Berezinzkii-Kosterlitz-Thouless transition in which the fermions acquire 
nonzero mass despite the absence of phase coherence. An analogous phenomenon is 
found in the physics of strong coupling or low carrier density superconductors. In these 
materials Cooper-pair formation occurs below a temperature T = T* whilst the U{1) 
symmetry remains manifest. The "local" gap modulus which is neutral under f/(l) ro- 
tations is nonzero, whilst the phase of the gap fluctuates violently this is known as 
the pseudogap phase. It is only below a temperature T = <^ T* that phase fluctu- 
ations cease and the material enters a true superconducting phase. This separation of 
the temperatures of pair formation and pair condensation has been known of for many 
years. It remains an open question, however, whether the dynamic generation of fermion 
masses without the breaking of chiral symmetry could be phenomenologically relevant 
in particle physics. 

In this paper we present a thermodynamic study of the 2 + Irf SU{2) (8> SU{2) 
Nambu— Jona-Lasinio (NJL) model on the lattice, with the aim of showing that these 
phenomena can occur in theories with more phenomenologically motivated symmetries. 

In the continuum, the model is described by the Euclidean Lagrangian density 

,2 



2Nf 



(1) 



where r = {ti,T2,T3) are the Pauli spin matrices which run over the internal SU{2) 
isospin symmetry and are normalised such that tr(raTa) = 26af3- The index i runs over 
Nf fermion flavors and is the coupling constant of the four-fermi interaction. 

The model is chirally symmetric under SU (2) SU (2) ijj — > {PlU +PrV)'iI), where 
U and V are independent global SU{2) rotations and the operators Pl,r = |(1 =t 75) 
project onto left and right handed spinors respectively. It is also invariant under U{l)v 
corresponding to a conserved baryon number. 

The theory becomes considerably easier to treat, both analytically and numerically, 
if we introduce scalar and pseudo-scalar fields denoted by a and vf respectively. The 
bosonised Lagrangian is 

C = + a + 2757f ■ f) + ^ (cr^ + 7f • 7f) 

- Nf + 

= +a + ij5^-T)iJi + -^tT<i>^<^, (2) 
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where the combination $ = cr + ivr ■ r is proportional to an element in the chiral group 
such that the model is invariant under the rotation $ V^U~^. 

Apart from any obvious numerical advantages this relatively simple model has various 
interesting properties: 

(i) the spectrum of excitations contains distinguishable baryons and mesons, i.e. the 
elementary fermions q and the composite qq states; 

(ii) at zero temperature and sufficiently strong coupling g'^ > g'^, chiral symmetry is 
spontaneously broken leading to a dynamically generated fermion mass equal to 
M(0) = (a) in the large-A^/ limit; the pion fields are the associated Goldstone 
bosons; 

(iii) for dimensionality 2 < d < 4 there is an interacting continuum limit at a critical 
coupling, which for d = 3 has a numerical value g'^/a ~ 1.0 to leading order in 
l/Nf with a lattice regularisation employed; 

(iv) the global symmetries are the same as those of 2 fiavor QCD, for which reason 
its four-dimensional version has a long history as an effective theory of QCD at 
intermediate energies. 

Having highlighted the model's zero temperature properties, let us now discuss what 
we know of its thermodynamics. To leading order in the effective potential has 

the same form as the discrete symmetry case with the replacement ^ + tt^. This 
implies that chiral symmetry remains broken up to a critical temperature = (M is 
the zero temperature fermion mass) and that symmetry restoration at this temperature 
is associated with the fermions becoming massless. This conclusion is expected to be 
valid only when Nj is strictly infinite, i.e. when the fiuctuations of the bosonic fields are 
neglected, since otherwise it runs foul of Coleman- Mermin- Wagner (CMW) theorem 0] 
which states that in 2 + 1 dimensional systems a continuous symmetry must be manifest 
for all T > 0. This behaviour may be understood by noting that, in the language of 
spin models, a domain wall separating regions of oppositely oriented magnetisation has 
an energy which does not increase with the size of the system. If a wall has length / 
and thickness t, the energy is ~ l/t, since the magnetisation is allowed to interpolate 
continuously [Q. This remains finite even as /,t oo. The physics is comparable, 
therefore, to that of the Ising chain, where a kink also costs only finite energy. At 
T = O"*" the magnetisation changes discontinuously, but there is no latent heat as the 
domains can grow to be very large. A simple estimate of the entropy reveals that the 
domain size grows exponentially with the inverse temperature. Studies of the critical 
properties of the 2 + Id NJL models with Z2 and U{1) 0, Q chiral symmetries 
at nonzero temperature have shown that the thermally induced phase transition of the 
Z2-symmetric model belongs to the universality class of the two-dimensional Ising model 
and that the U (l)-symmetric model has the same phase structure as the two-dimensional 
XY model. These results are in accordance with the dimensional reduction scenario, 
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which predicts that the long range behaviour at the chiral phase transition is that of 
the {d — 1) spin model with the same symmetry. This is because the IR region of the 
system is dominated by the zero Matsubara mode of the bosonic field. 

In the SU{2) ® S'f/(2)-symmetric model, one expects to observe fermion mass gen- 
eration in the absence of symmetry breaking at finite temperature, as in the case of the 
f/(l) model. The local amplitude p = y/a"^ + n ■ n is neutral under 0(4) rotations and 
may therefore be nonzero without breaking the symmetry. Hence, it is possible to have 
a dynamically generated fermion mass M(T) ~ p whose value may be comparable to the 
naive prediction of the large- iVj approach. This was demonstrated in d = (1-1-1) in ^j, 
in which results from an exactly soluble fermionic model were generalised to the ^7(1)- 
symmetric Gross- Neveu model. As was argued in 0, at T 7^ the fermionic spectral 
function is modified to a branch cut, implying that the propagating fermion constantly 
emits and absorbs massless scalars and hence has an indefinite chirality. The physical 
fermion may then be a superposition of positive and negative chirality and therefore 
be neutral under chiral transformations. It is in this way that one anticipates mass 
generation despite manifest symmetry. 

In this study we have carried out lattice simulations with Nf = 4. After briefly 
discussing the lattice formulation of the model, we present results of simulations at 
T = from which we have extracted the bulk critical exponents /3mag and u. For T 7^ 
we have attempted to study the issue of mass generation by measuring the so-called 
screening or "spatial" mass Ms{T) and the pole or "temporal" mass Mt{T). In order 
to better understand the relationship between the two we carried out a study of the 
lattice dispersion relation. Finally, simulation on lattices with various temporal extents 
helps to illustrate the nature of finite size and lattice discretisation effects. Measuring 
the pressure and making comparison with the Stefan-Boltzmann limit was found to 
elucidate the matter further. 

2 Lattice Formulation 

In its bosonised form, the model may be formulated on the lattice using the action 

^ — 2N 

-5 = E E [xiM[a, n]xi + C,M* [a, + — J2[a' + vf .vf) , (3) 

1=1 3; 9 X 

where C and ( each represent independent Grassmann-valued staggered fermionic 
fields defined on lattice sites x, and the auxiliary bosonic fields a and vf are defined on 
dual lattice sites x. The fermion kinetic operator is given by 



- ^yx-o) + E {6yx+u - Syx-u) 

i/=l,2 

+ bxyj:i^i^)s"'+H^Mi)-rn, (4) 

(x,x} 



4 



where {x, x) represents the set of 8 dual lattice sites neighbouring x, and the symbols 
riy{x) and e(x) are the phases (^— i)^o+---+2^i^-i g^^^j (^_]^~)a;o+xi+a;2 respectively, p, g are the 

internal SU{2) isospin indices, which we include explicitly. Making the replacement 
75 —>■ e{x), the symmetries defined in the previous section are still observed, i.e. the 
SU{2) ^SU{2) chiral symmetry remains exact |^, which is not the case in lattice QCD. 
With this replacement, the projection operators Pl,r —>■ Pe,o now project onto even and 
odd sub-lattices respectively. 

This choice of action corresponds to using a functional weig ht of det MtM. This 
doubles the fermionic degrees of freedom, but does allow one to perform simulations 
using a hybrid Monte Carlo algorithm, which has the advantage of being exact. Further 
doubling due to the use of staggered fermions in three dimensions leads to a total of 
Nf = 4N continuum fermion fiavors. Trajectory lengths were sampled from a Pois- 
son distribution with a mean of approximately 1.0. This ensures the ergodicity of the 
ensemble and reduces the autocorrelation of successive configurations. 

Simulations were carried out on lattices with x Lt sites, separated by lattice spacing 
a ~ 1/A, where A is the ultra-violet cut-off. In Euclidean spacetime, temperature is 
identified with the inverse temporal extent, i.e. T = l/{Lta). The temperature of the 
system may be altered, therefore, by varying either (3 = l/g"^, which amounts to varying 
the lattice spacing, or by varying Lt. The effective zero temperature limit is reached if 
one chooses Lt ^ Lg ^ I and cut-off dependent effects may be explored by varying j3. 
To reach the continuum limit one has to satisfy the following two conditions: A ^ T 
and A ^ m, where m is a mass scale inversely associated with the correlation length 
^ of the fluctuations of the order parameter. The condition A 3> T requires a lattice 
with sufficiently large Lj ^ 1. This condition not being met results in the sum over 
the Matsubara frequencies being truncated at too small a value. The second condition, 
A ^ m, implies that the coupling (3 should be tuned toward its critical value at which the 
correlation length ^ ^ 1. Furthermore, the condition Ls ^ C, should be satisfied in order 
to minimise finite size effects and obtain results which approximate the thermodynamic 
limit . 

3 Results 

3.1 Zero Temperature 

In order to study the chiral symmetry of the model first at T = and then at T > 
we work in the chiral limit, i.e. we choose not to introduce a bare fermion mass mo 
into eq. @. Without the benefit of this interaction, the direction of symmetry breaking 
changes over the course of the simulation such that S = ^ J2x '^(^) n = J2x ^{^) 
(where V is the lattice volume) average to zero over the ensemble. It is in this way that 
the absence of spontaneous chiral symmetry breaking on the finite lattice is enforced. 
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Instead, we define an effective order parameter 



^ = ( vs2 + n-n 



(5) 



which is a projection onto the direction of symmetry breaking for each configuration. 
Although this definition allows one to distinguish between the chirally broken and un- 
broken phases, it should be noted that the two definitions of the order parameter, $ 
and 



differ by a factor of l/yV as it is impossible to disentangle the chiral 
order parameter field in $ from the Goldstone modes. We also measure the fermion 
correlation function for temporal separation t; 



C{t) = At 



e 

C{t) 



Mt 



-M{Lt-t) 







if t 
if t 



odd, 



even, 



(6) 



which allows one to extract the fermion pole mass M. The fact that C{t = even) 
averages to zero over the ensemble reflects the fact that chiral symmetry is not truly 
broken on our finite size system, as manifest chiral symmetry implies that the Goo and 
Gee components of the staggered fermion propagator should vanish. 
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Figure 1: Effective order parameter and fermion mass as functions of j3 on 36^ and 48^ 
lattices. 

We simulated the model on a 36^ lattice with repeated simulations on a 48^ lattice 
near the critical coupling in order to detect and control finite size effects. Approxi- 
mately 1000 equilibrated trajectories were generated for each point. Fitting these results 
to the scaling forms $ = ai^P^"^^^ — f3Y'^''^ and M = 02 — PY we were able to ex- 
tract the critical exponents /5mag and u. As /5 ^ 0, the space-time lattice becomes large 
and coarse and the 0{a) lattice discretisation effects become significant. As /? — >• (3^^^^ 
the decreasing lattice spacing causes the lattice to become less coarse, but finite volume 
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effects become dominant. True scaling behaviour is limited, therefore, to a particular 
window, which we choose as 0.650 < j3 < 0.825 by demanding that results of the fits be 
stable with respect to the adding and subtracting of individual points. The measured 
values of the exponents were /J^ag = 1.12(2) and z/ = 1.10(2) with /J^"''^ = 0.961(5) and 
0.969(3) respectively. Data for $ and M versus P are plotted with the fitted scaling 
functions in Fig. |l]. 

Our results for z/ and /^mag are consistent with the 1/A^/ expansion results calculated 



where in our simulations Nj = 4. The fact that we observed next-to-leading order 
corrections, which were not observed in studies of the Z2 [§ and U{1) models, is 
because the 0{1/Nf) corrections in eq. (0) and eq. (|[) are slightly larger than in the 
models of lower symmetry. The fact that there are more light mesons that mediate the 
interaction in the S'[/(2) (S)S'[/(2)-symmetric model than in the Z2 and U{1) models may 
be the source for the larger deviation from the mean field values, which are /?mag = 1^=1 
for all three symmetries. More accurate values of the exponent could be extracted by 
performing detailed finite size scaling analyses in the vicinity of the critical point. 

In addition to the extraction of critical exponents, the studies presented here pro- 
vide a length scale, namely the fermion mass M, with which one can present nonzero 
temperature results in a cut-off independent fashion. This is performed in the next 
section. 

3.2 Nonzero Temperature 

In order to study the model at nonzero temperature we performed simulations on 
lattices with a constant temporal extent = 4 and varying spatial volumes Vs = 
30^50^100^150^ and 220^ The temperature was varied continuously by varying 
/3 < ~ 0.965, at which point a ^ and T ^ 00. 

The effective order parameter is plotted against temperature in Fig. ^. The tempera- 
ture is normalised by a factor of the zero temperature pole mass M(0), in order to make 
results cut-off independent as previously discussed. At high temperature, $ is consis- 
tent with zero whilst at low temperature there are significant finite size effects. These 
effects are a signal of the correlation length diverging as one approaches the transition 
at T = 0^ and the condition Ls ^ C, not being met. Once ^ ~ one is unable to 
observe a manifest symmetry. Our results are consistent with the expectation of chiral 
symmetry restoration in the infinite volume limit for T 7^ 0. 

This picture is supported by the results shown in Fig. |, in which the susceptibility 
of the order parameter 



m 



i 







(9) 
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Figure 2: Effective order parameter vs. normalised temperature for various spatial 
volumes. 



is plotted as a function of temperature. As spatial volume is increased, the pseudo- 
critical temperature Tc{Ls), i.e. the position of the peak in susceptibility, decreases, 
supporting a transition at T = O"*" in the infinite volume limit. Finite size effects are 
again very large in the low T regime, which is again consistent with the expectation 
that the correlation length diverges as T ^ 0. It is elucidating to compare our results 
for $ and x with those from simulations of the f/(l)-symmetric model 0, where it was 
shown that a Berezinskii-Kosterlitz-Thouless transition separates two chirally symmetric 
phases, one critical and one with finite correlation length. The main observations are: 
(a) in the U{1) model the finite size effects in the low T phase are larger than in the 
SU {2)^SU (2) model because the correlation length is infinite for any T < Tbkt] (b) the 
positions of the susceptibility peaks in the U{1) model converge to a nonzero T = Tbkt- 
In order to address the issue of mass generation with Lt = 4, we cannot directly 
extract the fermion mass from temporal correlators via eq. (^) as there are no longer 
enough temporal separations over which to carry out a fit. Instead, by studying the 
decay of spatial correlators, we extract the so-called screening or "spatial" mass Ms{T) 
via 



C{s) = if s = even. 



For free, massless fermions we expect this quantity to be the lowest fermion Mat- 
subara mode, given in the continuum by a;™"* = ttT, and on the lattice by wjf* = 
sinh~^ [sin(7r/L()]; cu]^^ ~ 0.658 for Lt = 4. In the low temperature phase of the large- iVj 
limit the "temporal" or pole fermion mass Mt{T) is related to the screening or "spatial" 
mass Ms(T) via = ul + M^, as the effects of the interaction are absorbed by the 
dynamically generated mass. If 1/Nf effects are taken into account, however, the rela- 
tion between and Mg can be less trivial. In particular, at nonzero temperature the 
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Figure 3: Susceptibility vs. normalised temperature for various spatial volumes. 

fermionic dispersion relation is expected to reflect the breaking of Lorentz invariance via 
a temperature dependent vacuum polarisation term n(/c,T): 



E'^{k, T) = M2(0) + + U{k, T). 



(11) 



Temperature effects can be absorbed into the thermal mass M{T) and a coefficient A{T) 
which can differ from one, implying a temperature dependent modification to the speed 
of light. The latter effect is expected to be significant for -C T, allowing one to write 



n(fc,T) = a'^{T)k'^ + 0{k'^). The dispersion relation becomes 

E\k,T) = M\T)+A\T)k\ 



(12) 



where = + 1. It is simple to show that the relationship between the screening and 
the pole mass is given by A^M"^ = + Uq, implying that a knowledge of Mg alone is 
not sufficient to reach a conclusion about the magnitude of the dynamical fermion mass 
Mf. This issue is investigated later in this section. 

Results from the measurement of the screening mass from the Lt = 4 simulations 
on various spatial volumes are presented in Fig. ^. The first observation to be made is 
that at low temperature Ms{T) is significantly larger than u]^^. It is an interesting aside 
to note that, as was the case with the U{1) model 0, the magnitude of the screening 
mass in our simulations is comparable to the values extracted from simulations of the 
Z2-symmetric model ^ where mass generation is ascribed to orthodox chiral symmetry 
breaking. Fig. ^ shows the mass to be size independent, even though, as we have seen 
from Fig. |^, the order parameter $ has large finite size effects for all T > 0. It is also clear 
from Fig. ^that Mg does not approach a value equal to tu^f* in the limit that T oo. We 
attempted to understand this by measuring the fermion mass on lattices with different 
temporal extents Lt = 8, 16 and 24 and for a range of values of the coupling (3. The 
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Figure 4: Screening mass vs. normalised temperature. The horizontal line represents 
the lowest Matsubara mode. 



screening mass measured on lattices with three different values of Lf are presented in 
Fig. 1^. The data are normalised by a factor of cujf'* to make it possible to plot the results 
on the same axes. At low temperature there seems to be no Lt dependence whilst at 
high temperature we see that the discrepancy between MgiT) and oj^q^ is reduced as Lt 
is increased. The effect of finite Lt is that only a finite number of Matsubara modes are 
present on the lattice, and when the fermion screening mass becomes comparable with 
the lowest Matsubara mode Ms{T) ~ cujf'*, discretisation artifacts become significant. It 
is clear that this effect is only apparent when comparison is made with large-Lt data. 

In order to gain insight into the effects of finite Lt and finite physical volume (which 
at nonzero temperature is translated into finite asymmetry ratio Lg/Lt) we studied the 
pole mass on different lattice sizes 8 x 48^, 16 x 48^, 16 x 96^ and 24 x 72^. The results of 
the pole mass normalised with respect to M(0) are presented in Fig. p. At this point we 
remind the reader that in order to increase T, we decrease the lattice spacing a, which 
is achieved by tuning the coupling (3 toward the bulk critical point Comparing 
curves with the same asymmetry ratios and different Lt, such as those extracted from 
data on lattice sizes 8 x 48^ with 16 x 96^ and 16 x 48^ with 24 x 72^, we see that the 
finite Lt effects become large as T increases. As in Fig. ^, this is due to the fermion 
screening mass becomes comparable with the lowest Matsubara mode. Similarly, the 
effects of finite spatial extent become more pronounced with temperature, since finite 
size effects become large when the fermion correlation length is comparable with the 
spatial extent of the lattice ^ ~ Mf^ ~ Lg. We have plotted the same data with a 
different normalisation in Fig. and fitted to the form Mt{T)/T = ai (T/M(0))"^. We 
chose to fit only to data with T/M(0)~0.5, as we know from Fig. ^ that in this region 
finite size and discretisation errors are not severe. This choice is justified by noting that 
in Fig. the data at large T deviate from the fitting line. We found a weak temperature 
dependence, with 0.68 < ai < 0.85 and —1.25 < 0^2 < —1-12 over the four fits. This is 
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Figure 5: Screening mass vs. normalised temperature. In order to compare different Lt 
we normalise by a factor of u?q^. 



consistent with — as T — > oo. In the future, it would be interesting to compare 
this result with next-to-leading order 1/Nf calculations. 

As previously discussed, in order to better understand the temperature dependence 
of eq. (|T2|) we have studied the fermion correlation function at nonzero momentum 
k = 27in/Ls {n = 0, 1,2, . . . ,Ls/4) in the spatial direction and used the energy E{k) 
extracted from eq. to map out the fermion dispersion relation. Due to the periodicity 
of the lattice in spatial dimensions and the doubling of fermion species, these dispersion 
relations are symmetric around k = 7r/2. We fitted data from the zero and nonzero 
temperature phases to the lattice free fermion dispersion relation O 



sinh^ E = sin^ k + sinh^ Mt. (13) 

Results from the hot phase on a 16 x 48^ lattice are shown with the fitted forms in 
Fig. ^. The masses extracted from fits to E{k) at both T = and T > are consistent 
with those extracted from the zero momentum correlators and yield values of A ~ 1 to 
within less than 3%. This implies that the principal physical effect of the hot medium 
is to generate a nonzero thermal mass, rather than to renormalise the speed of light. 
Furthermore, our results show that any extra factors not included in the free-field ansatz 
eq. ([T3|) are negligible. Similar results were extracted by fitting the dispersion relation to 
data from 16 x 96^ and 24 x 72^ lattices. The A ~ 1 result provided additional evidence 
that the nonzero screening mass extracted from our simulations with Lt = 4 implies that 
the fermion pole mass is nonzero. We see no modifications to the low momentum part 
of the dispersion relation, expected to have two branches at T > These branches, 



corresponding to two types of quasi-particle excitations (the fermion and the plasmino) 
are not visible in our simulations. 

Finally, to gain further insight into the thermodynamics of the model, we studied 
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the equation of state by measuring the pressure as a function of temperature using the 
integral method |TT|. For homogeneous systems the pressure P is given by 



P 



d\nZ 



dV 



InZ 



(14) 



where / is the free energy density. On the lattice the logarithm of the partition 
function is calculated from the expectation value of the bosonic part of the action 
5*608 = 2A^X]x[c"^(a;) + 7r(x) ■'7f(x)]. Since the derivative with respect to the bare coupling 
(5 is 

37 

-g^ = (S>. (15) 
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Figure 8: Fermionic dispersion relation vs. momentum k. The rescaled temperature 
T/M(0) decreases from 0.165 (top) to 0.735 (bottom). 

the physical free energy density can obtained from 



The free energy density is normalised by subtracting the vacuum contribution calculated 
at T = 0. The pressure was calculated from data on the 4x 100^ lattice whilst the vacuum 
expectation value of Sbos was calculated from data on 36^ and 48^ lattices. A lattice 
large- A^/ calculation of P, including the calculation of the lattice Stefan-Boltzmann limit, 
is outlined in Appendix A. This result calculated on a 4 x 100^ lattice is compared with 
the result extracted from simulations with Nf = 4 on the same lattice size in Fig. |^. 
The statistical error in the pressure extracted from simulation data is approximately 5% 
of the value of P for all T. It is striking that the value of P/T^ extracted from the 
simulations doesn't reach the Stefan-Boltzmann limit as T — oo, an effect found to be 
independent of the spatial volume. We believe that this is simply another manifestation 
of the effect that causes the screening mass not to approach Wq^* as T — oo. We interpret 
it as a renormalisation of discretisation artifacts by 1/A^/ corrections. Unfortunately, our 
data on lattices with > 4 were not sufficient to study the equation of state because the 
fact that as Lt becomes large the two contributions to eq. (^) converge, in conjunction 
with the normalisation factor of = Lij causes statistical errors to swamp any signal 
in the data. 

Drawing the information in this section together, our results provide clear evidence 
for the fermions acquiring nonzero mass in the absence of chiral symmetry breaking. In 
order to be certain whether mass generation switches off at some value of T = T* or 
whether it goes to zero asymptotically at T — oo, we would need to simulate the model 
on lattices with much larger temporal and spatial extents. 
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Figure 9: The pressure calculated from the lattice gap equation and from numerical 
simulations. In both cases the lattice size is 4 x 100^. 



4 Conclusions 

In this paper we have presented a thermodynamic lattice study of the 2 + 1 dimensional 
SU (2) (g) SU (2)-symmetric NJL model. We find the model to have manifest chiral 
symmetry for all T 7^ in accordance with the CMW theorem. As temperature is 
reduced, the susceptibility of the order parameter begins to diverge and is consistent 
with there being a phase transition at T = O"*". This should be contrasted with the 
predictions of the large- A'^y limit, in which fluctuations of the bosonic fields are neglected 
and chiral symmetry remains broken up to a critical temperature = ^^^^ limit, 

the breaking of chiral symmetry is associated with the generation of fermion mass in 
the standard way. We have demonstrated non-perturbatively, that in agreement with 
large- iVj there is mass generation at T 7^ in this model. 

By studying the so-called "screening" mass Ms{T), we have shown that at low but 
nonzero temperature Ms ^ ujq. With a temporal extent of Lt = 4, we observed that the 
spatial mass fails to approach ujq in the T 00 limit. This effect is also observed in the 
equation of state, where as T — >■ 00, the pressure fails to approach the Stefan-Boltzmann 
limit. This is due to the fact that when Lt is small there are only a small number of 
Matsubara modes available, an effect which becomes noticeable at high temperatures, 
where Ms ~ uq- Results from measurements of Ms on lattices with Lt > 4 have shown 
this to be a discretisation effect which is reduced when Lt ^ 1. By studying the fermion 
dispersion relation, we have also shown that any temperature dependent modification to 
the relationship between Mg and Mt is negligible, implying that we have unambiguous 
evidence that at T > we observe mass generation despite the lack of chiral symmetry 
breaking. It is interesting to note that the value of Ms{T) is of the same magnitude 
as it is in the model with Z2 chiral symmetry [^j, in which mass generation is ascribed 
to conventional symmetry breaking. To support our evidence of mass generation from 
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the measurement of Mg we have also measured the pole mass directly on lattices with 
Lt ^ 1. At high temperatures we observe large finite size effects as the correlation 
length of the low-mass fermion becomes comparable to the spatial extent of the system. 
At low temperatures, however, we observe that the mass moves smoothly toward its 
zero temperature value. As with the results from the U{1) model, we take this as 
further support for Witten's statement that when interpreted correctly, the large- 
Nf expansion is a reliable guide to the properties of the model. Although the CMW 
theorem forbids the large- A^j prediction of symmetry breaking at T > 0, the prediction 
of a dynamically generated fermion mass is adhered to. In this case the propagating 
fermion constantly emits and absorbs massless scalars and hence has indefinite chirality. 
As was argued in 0, the T 7^ fermionic spectral function is modified to a branch 
cut. In this model we observe mass generation without phase coherence, i.e. we have 
a pseudogap phase for T > 0. At this stage, however, we cannot say whether there is 
a phase transition at some temperature T* or the fermion mass goes asymptotically to 
zero as T — > 00. In order to study this regime (corresponding to a — >• 0) it is necessary 
to perform new simulations on lattices with larger Lt and Lg. 

To put these results into context, we note that the observation of mass generation 
without chiral symmetry breaking in 2 + Irf adds to the already rich and interesting 



phase diagram of the SU (2) ® SU{2) N JL model [0 . At high baryon chemical potential 



and zero temperature a complementary phenomenon is known to occur. Monte Carlo 
simulations have shown there to be a "thin film" superfiuid phase with long range phase 
coherence, but no mass gap is generated about the Fermi surface |]13| . Although similar 
in nature to the 2d XY spin model, a unique critical exponent seems to put the high-// 
phase in a new universality class. In 3 + 1 dimensions, the high-T low-/i phase shows 
conventional symmetry restoration via a T 7^ transition and the high-/i low-T phase 
appears to be that of a traditional BCS superfiuid. A study of the region at high-// 
and above the critical temperature, however, shows it to behave as if in a pseudogap 
phase [0]. Analogous to strong coupling superconductors there are two distinct critical 
temperatures: (a) T* below which there is manifest baryon number symmetry but no 
mass gap A, and (b) below which normal superfiuidity occurs. This is particularly 
interesting, as it provides evidence for mass gap generation without phase coherence in 
a phenomenologically relevant number of dimensions. 
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Appendix A: Large- A^j calculation of the Pressure 

In this appendix we outline a derivation for the pressure in the large-A^/ limit on the 
lattice. In section ^ we define our model in terms of staggered fermion fields which are 
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defined on sites x of tlie space-time lattice. In d dimensions, fermion doubling gives rise 
to 2'^ de grees of freedom per naive lattice flavor, which in the staggered formalism are 
interpreted as 2'^/ 4 "tastes"Q of Dirac spinor. This is easier to understand if we redefine 
the fields on a blocked lattice, with each site y associated with 2*^ sites of the original 



lattice |]T5|. The kinetic operator then tends to the continuum Dirac operator in the 
limit a ^ 0, with a pole in the propagator only at zero momentum. We apply this 
formalism, which is outlined for an odd number of dimensions in [Q, in the case that 
d = 3. Defining the 2 x 2 matrices Ta = r^'r^^r^^ and Ba = (-rl)^H-r2)^2(-r3)^^ 
we carry out a unitary transformation to new fields 

nriy) = j^.nriy) = TT'MA, y\ d^iy) = Ea bT^xU, y), (A-1) 

where y denotes a site on the "blocked" lattice with spacing 2a and A labels a site on a 
2^ cube with its origin at site y. Each site x of the original lattice corresponds, therefore, 
to a unique choice of y and A. The combination 



Qriv) = 1 ! , (A-2) 




may then be interpreted as a fermion field with spinor index a, "taste" index a and 
isospin index p. 

We identify the pressure with the negative of the free energy density — / = (T/Vs) In Z, 
where K is the spatial volume of the lattice. In the large-A^/ limit fluctuations in the 
bosonic fields are ignored and the partition function is given by 

Z = /dgdgexpf-f;Eg,Mt(fc)M(fc)g,')exp(-S'bos) 

= l^det M{k)j e-^^°% (A-3) 

where the fermionic action is written in terms of the fields q and q recast in Fourier 
space. The effective fermion kinetic operator as a function of fermion momentum is 

i 

^ i r , 
^(^) = E 9{(7M®l2®l2)sin2A:^ + (74®r;®l2)(l-cos2A;^)} 

+ (I4 ® I2 ® I2) S, (A-4) 

and the bosonic action is 22 

'S'bos = 5 • (A-5) 
9 



^We use the term "tastes" here to distinguish these 2 Kogut-Susskind flavors from our N lattice and 
Nf continuum flavors. 
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In the tensor products in eq. ( |A-4| ) the first matrix acts on spinor, the second on "taste" 



and the third on isospinor indices. S is the mean field equivalent of cr(x). Because the 
blocked lattice spacing is 2a, the Brillouin zone ranges from — 7r/2a to 7r/2a, meaning 
that the allowed momenta are 



= ^ n^ = 0,±l,...,±Lj2 for ^ = 1,2 
= ^^^^r % = 0,±1,...,±V2 for/i = 3, 



(A-6) 



where a has been set to 1. 

The determinant of eq. ([A-4[) is calculated as 



detM(fc) : 

so the unsubtracted free energy density is therefore 



1 ^ 

(sin2 2A;^ + 4sin^fcJ 



(A-7) 



1 —?N ?Ny'^ 

/ = -— - In Z = — - E In det M(A:) + (A-8) 

^s-^t ^s^t g 

As was done with the lattice data, we normalise this by subtracting the zero temperature 
contribution so that the pressure as a function of temperature is defined as 

P = -{fT-fo). (A-9) 

Finally, the Stefan-Boltzmann limit is found by calculating eq. ( |A-9| ) with S ^ 0, which 
corresponds to T ^ oo. 
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